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^ \ Systematic experiments on stress-induced phase transitions in thin SMA 

^ ' structures in literature have revealed two interesting instability phenom- 

■ - - ' ena: the coalescence of two martensite-austenite fronts leads to a sud- 

den stress drop and that of two austenite-martensite fronts leads to a 
sudden stress jump. In order to get an insight into these two phenom- 
ena, in this work we carry out an analytical study on the stress-induced 
phase transitions in a thin SMA layer (a simple structure in which the two 
coalescence processes can happen). We derive a quasi-2D model with a 
non-convex effective strain energy function while taking into account the 
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rate-independent dissipation effect. By using a coupled series-asymptotic 
expansion method, we manage to express the total energy dissipation in 
terms of the leading-order term of the axial strain. The equilibrium equa- 
tions are obtained by maximizing the total energy dissipation, which are 
then solved analytically under suitable boundary conditions. The analyti- 
cal results reveal that the mechanism for such instabilities is the presence 
of "limited points", which cause the switch of nontrivial solution modes to 
trivial solution modes. Descriptions for the whole coalescence processes of 
two fronts are also provided based on the analytical solutions, which also 
capture the morphology varies of the specimen. It is also revealed the key 
role played by the thickness-length ratio on these instabilities: the zero 
limit of which can lead to the smooth switch of nontrivial modes to trivial 
modes with no stress drop or stress jump. 

1 Introduction 

Systematic experiments on stress-induced phase transitions in thin SMA 
structures, such as wires, strips and tubes, have been carried out (Lexcellent 
& Tobushi 1995; Shaw & Kyriakides 1995, 1997; Sun et al 2000; Favier et al 
2001; Li & Sun 2002; Feng & Sun 2006). In general, the whole experimental 
process can be divided into three stages. In the first stage, the deformation of 
the specimen is the elastic deformation of initial phase. In the second stage, the 
product phase first nucleates at some special site, and then propagates gradually 
along the specimen. Finally, the phase fronts coalescence and the whole length 
of the specimen has transformed to the product phase. In the third stage, the 
deformation of the specimen is the elastic deformation of the product phase. It 
was also found that the measured engineering stress-strain curves have some key 
features, e.g., the nucleation stress peak (for the loading case) and stress valley 
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(for the unloading case), the stress plateaus, the rate-independent hysteresis loop 
and so on. 

In this paper, we shall focus on the instability phenomena induced by the 
coalescence of two phase fronts. During the experiments, these instability phe- 
nomena can be easily identified from both the stress-strain curves and the surface 
morphology of the specimens (Shaw & Kyriakides 1997; Feng & Sun 2006). Sys- 
tematic experimental results have obviously shown that the coalescence process is 
inevitably accompanied the varies of the stress value and the surface morphology 
of the specimen. In general, the coalescence process has the following procedure 
as described in Feng & Sun (2006). In the loading (unloading) the total 

elongation increases (decreases), two martensite-austenite (austenite-martensite) 
phase fronts move towards each other. The austenite (martensite) region sand- 
wiched between the phase fronts becomes more and more narrow. Eventually 
these two fronts start getting touch each other and the coalescence process takes 
place. It is clear the coalescence process is a dynamic process and seemingly 
uncontrollable in the sense that it occurs very rapidly (this may indicate that an 
instability occurs). During the coalescence process, the local configuration of the 
specimen transforms from an inhomogeneous mode to a homogeneous mode and 
the corresponding stress value has a rapid drop for the loading case and a rapid 
jump for the unloading case. 

In order to get an insight into the phase fronts coalescence phenomena, in 
this work we carry out an analytical study on the stress-induced phase transitions 
in a thin SMA layer (a simple structure in which the two coalescence processes 
can happen). We shall derive a quasi-2D continuum model with a non-convex 
effective strain energy function while taking into account the rate-independent 
dissipation effect. The starting point of our model is the formulation of Rajagopal 
& Srinivasa (1999, 2004) (also see the formulation of Sun & Hwang 1993a, b). 
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In fact, we have already proposed a quasi-3D continuum model earlier to study 
the phase transitions induced by extension in a slender SMA cylinder (see Wang 
& Dai 2009), which has almost the same formulation with this current model. 
Thus, we shall only give a brief introduction to the derivation procedure in this 
paper and the detailed derivation procedure can be found in Wang & Dai (2009). 

In the experiments, typically the thickness-length (or width-length or radius- 
length) ratio of the specimen is of order O(10~^). As a result, one might think 
that the deformation along the lateral direction can be neglected and the layer 
can be treated as a one-dimensional object. However, sometimes a purely one- 
dimensional model appears to be not sophisticated enough to capture some key 
features observed in experiments. Some explanations have already been given 
in Dai & Cai (2006). In this paper, the total elastic potential energy of the 
layer will be considered based on a two-dimensional setting. Starting from the 
two-dimensional governing system and by using the coupled series-asymptotic 
expansion method (Dai & Cai 2006; Cai & Dai 2006), we manage to express the 
total elastic potential energy of the layer in terms of the leading order term of 
the axial strain (cf. (3.41)). Although the final expression of the total elastic 
potential energy is one-dimensional, it takes into account the lateral deformation 
and has the higher-dimensional effects built in. 

To describe the hysteretic behavior during the phase transition process, one 
also needs to consider the inevitable mechanical dissipation effect. In this paper, 
the mechanical dissipation effect will be considered in a purely one-dimensional 
setting, i.e., we neglect the influence of the radial deformation on the mechanical 
dissipation. A specific constitutive form of the rate of dissipation function is 
adopted in our model (cf. (2.6)). Based on the phase transition criteria and the 
criterion of maximum rate of dissipation, we propose the evolution laws of the 
phase state variable for the purely loading and purely unloading processes. After 
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doing some further analysis, we derive the one-dimensional expressions for the 
mechanical dissipation functions in terms of the axial strain (cf. (4. 5)- (4. 8)). 

With the expressions of the total elastic potential energy and the mechanical 
dissipation functions, the equilibrium configurations of the layer for the purely 
loading and purely unloading processes can be determined by using the principle 
of maximizing the total energy dissipation. By using the variational method, 
we derive the equilibrium equations, which are then solved analytically under 
suitable boundary conditions. It will be seen that the solutions obtained show 
qualitatively agreements with the experimental results. 

Based on the analytical solutions obtained and by using the limit-point in- 
stability criterion, we further consider the phase fronts coalescence process. It 
is revealed that during the coalescence process, the configurations of the layer 
switched from the nontrivial solution modes to the trivial solution modes, which 
is caused by the presence of the "limit points". The morphology varies of the 
layer and the accompanying stress drop/jump during the coalescence process can 
be described. The influence of the thickness-length ratio of the specimen on the 
coalescence process is also studied. It will be shown that the zero limit of the 
thickness-length ratio can lead to the smooth switch of nontrivial solutions to 
trivial solutions with no stress drop or stress jump. 

This paper is arranged as follows. In section 2, we give a simple introduction 
to the formulation of Rajagopal & Srinivasa (1999, 2004) and propose the evolu- 
tion laws of the phase state variable for the purely loading and purely unloading 
processes. In section 3, we formulate the field equations by treating the thin layer 
as a two-dimensional object. By using the coupled series-asymptotic expansion 
method, we express the total elastic potential energy of the layer in terms of the 
asymptotic axial strain. In section 4, we study the mechanical dissipation effect. 
After some analysis, we express the mechanical dissipation functions in terms of 
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the axial strain. In section 5, we derive the equihbrium equation by using the 
principle of maximizing the total energy dissipation. Then, we construct the an- 
alytical solutions for an illustrative example. In section 6, we further study the 
instability phenomena induced by phase fronts coalescence. We try to give some 
descriptions and explanations for the origin of the instability during the coales- 
cence process, the accompanying stress drop/jump and the morphology varies of 
the specimen. We also consider the size-effect of the specimens on the coalescence 
process. Finally, some conclusions are drawn. 

2 Preliminaries 

In this section, we shall give a simple introduction to the formulation of 
Rajagopal & Srinivasa (1999, 2004), which is the starting point of our present 
model. 

First, based on the balance laws and the local form of the entropy production 
equation, one can obtain the reduced energy-rate of dissipation relation as 

-(^ + Tr/) + tr(SF) =C. (2.1) 

In equation (2.1), \1/ is the Helmholtz free energy per unit referential volume, T 
is the absolute temperature, rj is the entropy per unit volume in the reference 
configuration and ( is the rate of mechanical dissipation, F is the deformation 
gradient tensor and S is the nominal stress tensor. The superposed dot indicates 
the material time derivative. 

To describe the phase transition process, one also needs to adopt the phase 
state variable a, which represents the volume fraction of the martensite phase in 
the reference configuration. It is assumed that S and r] depend on the variables 
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{F, a, T}, while the rate of mechanical dissipation ( depends on {T, a, a}. With 
these constitutive assumptions, the reduced energy-rate of dissipation relation 
(2.1) can be rewritten as 



tr((E - ^)F) - + r7)T - = C, (2.2) 



Equation (2.2) must hold for all F and T, thus one can arrive at the following 
constitutive equations: 

_ c^^(F,a,T) 
^ - dF ' ^^-^^ 

' dT ' ^ ^ 



and 



9^fF a T) 

^ ' '^-d = C. (2.5) 



da 

Once the specific forms of the functions ^ and ( are given, the nominal stress 
tensor S and entropy rj can be calculated from (2.3) and (2.4). Equation (2.5) 
provides the equation for the evolution of phase state variable a. As we just 
consider the isothermal responses of SMA materials in this paper, the temperature 
T will be considered as a given constant in the sequel. 

In the paper of Rajagopal & Srinivasa (1999), the following constitutive form 
of the mechanical dissipation rate ( was proposed 



A+(a)|d|, if d > 0, 

0, if d = 0, (2.6) 

A~(a)|d|, if d < 0, 



where A'^{a) > and A [a) > referred to as the forward and backward 
dissipative resistances, respectively. An important consequence of (2.6) is that 
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the response of the material is rate-independent, i.e., the stress-strain curve is 
the same irrespective of the speed with which the test is conducted. 

Substituting (2.6) into (2.5) and using the fact that the resulting equation 
must be valid for all admissible processes, one can arrive at the following phase 
transition criteria 

-A-(a)<-— <A+{a)^a = 0, (2.7) 
oa 

I - A (a), if a < 0. 

Condition (2.7) represents the fact that as long as —d'^/da lies between —A~ 
and A'^, no phase transition can take place. On the other hand, whenever a ^ 0, 
—d'^/da must equal to one of the dissipative resistances. At the two critical 
points —d'^/da = A'^ or —d'^/da = —A~, the response of material may not 
be unique (Rajagopal & Srinivasa 1999, 2004). To fully determine the com- 
mencement and cessation of the phase transitions, one need to use the criterion 
of maximum rate of dissipation (Rajagopal & Srinivasa 1998, 1999). Intuitively 
speaking, the maximum rate of dissipation criterion states that if the material is 
capable of responding in many different modes with different rates of dissipations 
(including a mode that is non-dissipative), then the actual mode will be the one 
with the maximum rate of dissipation. 

In general, there does not exist a one-to-one mapping between F and a. 
Based on the phase transition criteria and the criterion of maximum rate of 
dissipation, we find that for a given deformation gradient F, there may exist 
two phase state values and such that — 9\1/(F, ap)/(9a = A'^{ap) and 
— 9\1/(F, a;p)/9a = —A~{a^), then the actual phase state a should satisfy < 
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In this paper, we shall only consider two kinds of loading patterns: the 
purely loading process with the initial state of the cylinder composed of austenite 
phase and the purely unloading process with the initial state of the cylinder 
composed of martensite phase. In the figure of the stress-strain response, the 
curves corresponding to these two kinds of loading patterns form the entire 'outer 
loop'. 

From the experiments, it was found that the deformation process at any 
material point can be divided into the elastic deformation process of the austenite 
phase, the elastic deformation process of the martensite phase and the phase 
transition process. The experimental results also reveal that the profile of the 
phase transition (localization) region has a stable form during the phase transition 
process. Thus, it is reasonable to impose the following assumption: 

Assumption 1. In a quasi-static purely loading or purely unloading process, 
for the material points located in the phase transition region, it is satisfied that 
q; > (for the loading process) or d; < (for the unloading process). 

With this assumption, we can propose the following relationship between a 
and F for the purely loading and unloading processes, respectively: 



if - 



9^(F,a) 



< A+(0) 



da 



Q = 



a+(F) = <^ 1, 



if 



di>{¥,a) 



>A+(1), 



da 



a=l 



a 



if 3 a+ G (0,1), s.t. 



a^(F, a) 



A+ia+). 



da 



(2.9) 
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a-(F) = < 



0, 

1, 

a' 



a^(F,a) 

da 
9^(F,a) 



|a=o < -^"(0), 



F' 



da 

if 3 ttp G (0, 1), s.t., 



1 > -^-(1), 

a^(F,a) 



da 



-A (ap). 
(2.10) 



3 The total elastic potential energy 

In this section, we shall derive an "effective" one-dimensional expression for 
the total elastic potential energy of the layer, which takes into account the higher- 
dimensional effects. 

We consider the symmetric deformation of a thin SMA layer subject to a 
static axial force at two ends. This layer can be considered as a cross-section cut- 
ting along the axial direction of a thin SMA strip or tube. Although here we just 
consider a very simple model as we treat this layer as a two-dimensional object, 
it will be seen that this simple model can capture some important experimental 
features. 

Assume that in the stress-free configuration, the layer has thickness 2a and 
length I, where a/l = 6 <^ 1. We denote {X,Y) and {x,y) the coordinates of 
a material point of the layer in the reference and current configurations, respec- 
tively. The finite displacements can be written as 

U{X, Y) = x(X, Y) - X, W{X, Y) = y{X, Y) - Y. (3.1) 

As here we only consider the symmetric deformation of this layer, thus U and V 
have the following properties 



U{X,-Y) = U{X,Y), W{X,-Y) = -W{X,Y). 
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Then the deformation gradient tensor F is given by 



F = (1 + Ux)ei ® El + Uyei ® E2 + Wxe2 ® Ei + (1 + WY)e2 ® E2, (3.2) 

where {Ej}j=i_2 and {ej}j=i^2 are the orthonormal base vectors in the reference 
and current configurations. 

Based on (2.3), (2.9) and (2.10), we obtain that the nominal stress tensor T, 
for the purely loading or purely unloading process should be given by 



^ - gp \a=a±{F)--^ (*J, (3.3) 



where are tensor-valued functions. It will be assumed that there exist some 
"effective" strain energy functions \E'^(F) such that 

. G-,F) . (3.4) 

In this section, the purpose is to derive an "effective" one-dimensional ex- 
pression for the elastic potential energy of the layer. Thus, we only consider the 
nondissipative case in this section, i.e., we assume A±(a) = 0. From (2.9) and 
(2.10), it is easy to see that in the nondissipative case, the phase state functions 
satisfy a+(F) = a^(F) =: a{F). Further by using (2.9), (2.10) and (3.3), we 
have 

5^(F) = ||(F,a)U=„(p) 

d4f^^ da(F) /o p-n 

_ d^{F,a{F)) 
~ dF ■ 

The final equation in (3.5) is due to the fact that if F takes value in the phase 
transition domain, we have —^{F,a)\a=a{F) = 0, otherwise a{F) should be a 
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constant function of F. 

From (3.5), we can see that the "effective" strain energy function in the 
nondissipative case should be given by 



*(F) = ^(F,a(F)). (3.6) 

In the sequel, we refer \E^(F) as the elastic potential energy function. To go 
further, we propose another important assumption: 

Assumption 2. In the nondissipative case, SMA material can be considered 
as some kind of isotropic hyperelastic material. In this case, the elastic potential 
energy \E'(F) only depends on the two principle stretches Ai, A2 of F; that is 
* = ^(Ai, A2). 

If the components of F — I are relatively small, it is possible to expand the 
nominal stress components in term of the strains up to any order. The formula 
containing terms up to the third order material nonlinearity is (cf. Fu & Ogden 
1999) 



~ 0']ilk'^kl + r)^%iknm'^klK-rnn + r^O.^lknmqp'^kll^mnl^pq + 0{\Kst\^), (3.7) 



where the components of the tensor F — I and 



a%knmqp dF^^QF^ldEmndF^, '^"^ 



are incremental elastic moduli, which can be calculated once a specific form of 
elastic potential energy is given. From the formula (3.7), we can obtain the 
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nominal stress components E^j. For example, 



hi =T2Wy + TlUx 

I (3i 

+ -iw^iUlWY + e2Wl + w^ulUx + ?>e^wlUx + 3^2wvf/l 

6 

+ diUl + QOjUyWyUx + QOJJyUxWx + 3e6W^yW^l + 3e3f/xW"l), 



where r,, ?7j and 0^ are some elastic moduli, whose formulas are given in appendix 
A. Owing to the complexity of calculations, we shall only work up to the third- 
order material nonlinearity. In the experiments, the maximum strain is less than 
10%, and such an approximation is accurate enough, certainly not worse than 
the trilinear approximation adopted in many works. 

For a static equilibrium configuration, the nominal stress tensor S satisfies 
the field equations 

Div(S) = 0, 
which yields the following two equations 



dX dY 



0, (3.9) 
0. (3.10) 



dX dY 

We consider the traction-free boundary conditions at y = ±a, i.e. 



^21 1 y=±a 



0, S22|y=±a = 0. (3.11) 



Equations (3.9) and (3.10) together with (3.11) provide the governing equations 
for two unknowns U and W . 

We can also expand the elastic potential energy \l/ in terms of the strains up 
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to any order. The formula containing terms up to the fourth-order nonhnearity 
is given by 

* =linU^ + nW^ + 2t2WyUx + nU^ + 2t3UyWx + TsW^) 



+ r/if/|. + eriiUYWyWx + eri^UyUxWx + Sfj^WyW^ + Sri^UxW^) 
+ ^(^sf/y + GOsU^W^ + 9iW^ + m^U^WyUx + Ae2W^Ux 
+ mUyUl + Gd^W^Ul + 4^2H/'y t/i + ^if/x + 4^9f/y VTx 

+ 12^4[/yiy2p^x + 2Ae^UyWyUxWx + 12^4f/y f/^-l^x + QOiqUIWI 

+ 6^3iy^iy| + i2e6W^yt^xW^l + Qe^u^wl + 4e9f/yiy| + e^w^). 

(3.12) 

It can be seen that (3.12) has a very complex form. To obtain the total elastic 
potential energy, one need to calculate the integration of (3.12) over the whole 
layer, which will become more complex. 

Here, we shall adopt a novel approach involving coupled series-asymptotic 
expansions to obtain an asymptotic expression of the total elastic potential energy 
of the layer. A similar methodology has been developed to study nonlinear waves 
and phase transitions in incompressible materials (see Dai & Huo 2002, Dai & 
Fan 2004, Dai & Cai 2006, Cai & Dai 2006). 

Based on the symmetry of the problem (see the equation below (3.1)), we 
first introduce the important transformation 

W{XX) = Yw{XX). s = Y^-. (3.13) 
From the properties of U and W , it is easy to see that U and w can be considered 
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as functions of the variable s. Then we adopt the scahngs: 



s = I'^s, X = Ix, U = hu, w = —w, e = y. 




(3.14) 



where /i is a characteristic axial displacement, and e (equivalent to a small engi- 
neering strain) and u (square of the half thickness-length ratio) are regarded to 
be two small parameters. 

Substituting (3.13) and (3.14) into (3.9) and (3.10), we obtain 

2T3Us + {t2 + T3)W^ + TiUxx + S^ATsUss + (2r2 + 2t3)W:cs) H = 0, (3.15) 



Here and hereafter, we have dropped the tilde for convenience. The full forms of 
(3.15) and (3.16) are very lengthy and can be found in appendix B. Here we just 
present the first free terms. Substituting (3.13) and (3.14) into the traction- free 
boundary conditions (3.11), we obtain 

2t3Us + T^Wx + e[2r]2,wus + ^rjsUsUx + v^ww^ + V^UxW^ 

+ ]^9iw'^Wx + djwu^Wx + ]^OAul.Wa: + s^{A9-iUswl + 2diwlwx) (3.17) 
4 

+ s{-9sul + AO^wUgWs + AOqUsWsUx + 29quIwx + 264WWsWx 
3 

+ 2e-jWsUxW^ + OiqUswI + \9<^wl)]\s=u = 0, 



GTiWs + {2t2 + 2t3)Uxs + TsWr^r, + sAtiWss H = 0. 



(3.16) 
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TiW + T2Ur, + s2tiWs + e[^T]iw'^ + S^2r]iwl + 7]2WU^ + ^?72M^ 

+ s{2r]3ul + 2r]iWWs + 2r]2WsU^ + 2r]4UsW^ + ^r^s^^^)] 

D O 2 Z 

+ s^(46'3M^w^ + 26*1 w;u7^ + 26*2 + 46'4U,M;,ty^. + O^Wswl) 



(3.18) 



+ s{26swu1 + Oivo^Ws + 26quIux + 292WWsUx + 65Wsul+ 

2d^WUsW^ + 2djUsU^W^ + ^6^3^^^ + ^6'6MxW^)]|s=i/ = 0. 

In the nonlinear PDE system (3.15)-(3.18), the two unknowns m and de- 
pend on the variable x, the small variable s and the small parameters e and i^. 
As long as we assume that u and w are smooth enough in s, we can expand them 
in terms of the small variable s: 

u{x, s] e, v) = Uo{x; e, u) + sUi{x; e, u) + s'^U2{x; e,u) -\ , (3.19) 

w{x, s; e, u) = Wo{x; e, u) + sWi{x; e, u) + 5^1^2(2;; e, i/) H . (3.20) 

Substituting (3.19) and (3.20) into the traction-free boundary conditions (3.17) 
and (3.18), then by neglecting the terms higher than 0(ei^, e^), we obtain 



2r3f/i + tsWq^ + e{2r]3UiWo + 2r]3UiUo^ + v^WoWq^ + r]4UoxW^ 



Ox 



+ \9,W^Wo. + 9jW,UoxW,, + \94Ul,W,,) + ez/(4r/3f/2W^o (3-21) 
+ Qti^UiWi + 4r/3t/2f/ox + 2r/3f/if/ix + 3r/4l^il^ox + r/4f/ixW^oa; 
+ r?4W^oW"ix + TliUo^Wi^) = 0, 
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tiWq + T2U0X + e{]^T]iW^ + Ti2WoUox + \v2Ul^) + //(SriVTi 

6 2 2 6 (322) 

+ 2r/4t/il^ox + ^^3W^ox) =0. 

It can be seen that (3.21) and (3.22) involve the five unknowns: t/o, Wq, Ui, Wi 
and 1/2- Thus, to form a closed-system, we need to find another three equations 
which contain and only contain these five unknowns. 

Substituting (3.19) and (3.20) into (3.15), the left hand side becomes a series 
in s and all the coefficients of s" (r;,=0, 1, 2, ■ ■ ■ ) should vanish. Equating the 
coefficients of s° and to be zero yield that 

2r3f/i + (ra + r3)l^ox + nUo^^ + e{27]sUiWo + 
+ (?72 + rii)WoWQ^ + (?72 + ri4)UQ^WQ^ + ?72W^ot^oxx + ViUoxUqxx) 
+ e\esUiW^ + 2e6f/il^of/o. + OsUiUl + (y + y)Vro'W^ox (3.23) 

+ (^5 + Oj)WoUo,Wo, + (| + j)u^,Wo. + \e,w^u,,, 

12r3f/2 + (3r2 + 3x3)1^1, + nU,^,^ + e{12r,-,U2Wo + 18r/3^7il^i 

+ 12r/3t/2f/o. + 10r]sU,Ui, + (3r/2 + 9r],)W,Wo, + (r/2 + 5r],)Ui,Wo, 

+ (3r/2 + 3rj^)WoWi, + (3r/2 + 37/4)^/0,1^1, + 37]2W,Uo,, + ViUi.Uq,, 

+ + r/if/o,?7i,, + 2r/4?7iiyoxx + VsWo.Wo,.) + e'(4^8?7i^ (3-24) 

+ 6esU2Wo^ + 1803^/1^^01^1 + 129qU2WoUo, + meU^WiUo^ 

+ 603t/2t/o,2 + m,U,WoU,, + m,U,Uo,U,, + {20, + 609)f/i2W^ox 

+ (302 + 904)1^0^^1 H^ox + (305 + 9e7)WiUo.Wo, + (05 + 507)l^ot/ixW^o. 
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+ {02 + 5e^)Uo.UiM. + (3^10 + 2ey)UiWo^^ + + 



1 2 1 



+ 92WoUo,Ui,, + ^^it/ox't/i,., + 2erUiWoWo.. + 29^UiUo,Wo., 

Similarly, substituting (3.19) and (3.20) into (3.16) and equating the coefficient 
of to be zero yields that 



GnWi + (2r2 + 2r3)f/i. + nWo,, + e(4r/3f/i' + 6riiWoWi 
+ {2r]2 + 2r]^)WoUi, + (2r/2 + 27]^)Uo,Ui, + 6r]^UiWo, + 2r/3Vrox 
+ 2r/4t/if/o^^ + r/3Wo:,f/oxx + Vs^oWq^^ + r/3t^oxW^o,;,; + 6r/2Vrif/o^) 

+ {62 + ^4)l^o't/lx + (2^5 + 2e7)iyot/oxf/l. + (^2 + 0^)UlUi, 



(3.25) 



+ 6^4f/iW^oW^ox + QOrU^Uo^Woa: + 2^3W^oW^ol + 2eeUoa:Wl 

+ 2ejUiWoUo,, + 29JJiUMxx + OqWoWMxx + ^3^^0xVro,?7o,a= 
+ \e^Wo'Woxx + ^eW^of/oxP^O.. + ^^3f/ox^Oxx) = 0. 

Now the nonlinear PDE system (3.15)-(3.18) has been changed into a one- 
dimensional system of differential equations (3.21)-(3.25) for the unknowns f/o, 
Wq, Ui, Wi and U2. Based on (3.23)-(3.25) and by using a regular perturbation 
method, we can express Ui, Wi and U2 in terms of Uq and Wq. The results are 
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nven below: 



Ui ={-\ - 7^)Wox - ^Uoxx + e^aiWoWox + ci2UoxWox 

+ a^W^U^xx + a^UoxUoxx) + e^a^W^Wox + aeWoUoxWo^ ^^'^^^ 
+ ayt/o^Wox + asW^Uoxx + adWoUoxUa^x + aio^oo^^oxx), 



^1 =(^7^ + TT^Woxx + + ^Woxxx 



A. 



+ e(aiiiyol + auWoxUoxx + aisf/ox^ + auWoWoxx 

+ Ol5^0xWo2:x + CLiqWqUoxxx + ^^17f^x•^^0a;a:a:) , 
1 T*^ T 

+ ( 1 ^ ^ )Uoxxxx + efflisW^oxW^oxi 

^ 24 24r| 24r| 12^ ^ (3.28) 

+ CtigUoxxWoxx + 0.2oWoxUoxxx + a2lf^0xxf^0xa:a: + <^22WqWqxxx 
+ a23f^0xW^0xxx + 0'24WoUqxxxx + 0,25UoxUqxxxx) , 

where (z = 1,2, ■ ■ ■ ,25) are material constants related to the elastic moduli, 
whose expressions can be found in appendix C. Substituting Ui, Wi and U2 into 
(3.21) and (3.22) and omitting the higher order terms yield the following two 
equations with only two unknowns Uq and Wq: 

- r2Wox - tiUqxx + ^{-m^oWox - 'n2UoxWox - 'n2WoUoxx - ViUqxUqxx) 

+ '^{^^^T^Woxxx + '^^^^Uoxxxx) + e\-\e2W^Wox ' O^W^UoxW^x 
2 

- \d2UlxW0x - \e^W^Uoxx - e2WoUoxUoxx - lo^ulUoxx) (3-29) 
+ eu{biWoxWoxx - hUoxWox - hWoxUoxxx + hUoxxUoxxx - hWoWoxxx 

+ beUoxWoxxx — brWoUoxxxx — bgU^xUaxxxx) = 0, 
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^ 2 2 6 2 2 ^330) 

D 

^ bisUoxWoxx — bii^oUoxxx — bi^UoxUoxxx) = 0. 
Integrating (3.29) once, we obtain 

C - T,Wo - nUox + ei-^V2W^ - V2W0U0X - ^ViUl) + i^i^^^^Woxx 

+ ^^^^Uoxxx) + e'i-h.wi^ - h.W^Uox - \e2W,ul - U,ul) 

+ eu{-bi6WQ^ - bnWoxUoxx - bisU^^^ - biQWoWoxx + b2oUoxWoxx 

— b2lWoUoxxx — b22UoxUoxxx) = 0, 

(3.31) 

where C is the integration constant. The coefficients bi {i = 1,2, ■■■ , 22) in (3.29)- 
(3.31) are also some constants related to the elastic moduli, whose expressions 
can be found in appendix D. 

It is important to find the physical meaning of C, since to capture the in- 
stability phenomena observed in the experiments, one needs to study the global 
bifurcation as the physical parameters vary. For that purpose, we consider the 
resultant force P acting on the material cross-section that is planar and perpen- 
dicular to the X-axis in the reference configuration, and the formula is 




(3.32) 



The formula for the nominal stress component En has already been given in (3.8). 
Through the transformations (3.13)-(3.14) and the series expansions (3.19)-(3.20), 
we can express En in terms of Uq, Wq, Ui, Wi and f/2. Further by using (3.26)- 
(3.28), we get an expression of En which only depends on Uq and Wq- Then, 
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carrying out the integration in (3.32), we obtain that 



1 1 

P =2ae[r2Wo + nf/o^ + e{-r]2W^ + r/aW^o^^Ox + ^ViUox) + 

D / 



I 



(3.33) 



+ hi^WoWQxx - hoUoxWoxx + hlWoUoxxx + b22UoxUoxxx)]- 

Comparing (3.31) and (3.33), we have C = P/{2ae). 

By using the two equations (3.29) and (3.30), we can express Wq in terms of 
Uq. First, from (3.29), we have 

Wqx = —{-nUoxx + e{-r]2WoWox - r]2UoxWox - V2W0U0XX - ViUoxUoxx))- 

T2 

(3.34) 

Substituting (3.34) into the term i'{^T2Woxx + ^TiUqxxx) of (3.30), we obtain the 
following equation 



I z o 



+ \mlu,x + \m,ul + -e2Ulx) + ev^-^f^^{-Wl - (3.35) 

2 z Z[ti — T2) 

— 2WqxUoxx — ^O^Oxx ~ UqxWqxx " ^qUoxxx ~ UqxUoxxx) = 0. 

From (3.35) and by using a regular perturbation method, we obtain 



Wo = - —Uox + eaiUl^ + e^a2Ul^ 
1 1 

+ - ^)UL + (2 - ^WoxUoxxx), 



(3.36) 
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where 



T"2 '7i + ^1 ^72 - 2rir2?72 



2rl 



1 

+ 3riV2(-?7i?72 - 2?7^ + ri6'5)). 



By using the coupled series-asymptotic expansion method introduced above, 
we can derive a one-dimensional asymptotic expression for the total elastic po- 
tential energy of the layer. The elastic potential energy per unit referential length 
of the layer is 

^> = ! mdY. (3.37) 



The formula for the elastic potential energy \1/ has already been given in (3.12). 
Through the manipulations we just described above, we obtain 



2 6 3(ri - ra) 3(ri - T2) 
K^Lx + 7rT-2VroVroxx + 7T(-ri - 2r2)f/o.iyoxa: 



3(t-i - r2) 6 6 

11 1 

+ -tiWqUq^^^ + g(-2ri - r2)f/oxf/oxxx) + e^(^^iW^o (3 33) 



CgUo^Woxx + CioWqUqxxx + Ci i f^Ox ^Oxxx + CuUq^Uoxxx)] 



where Cj (z = 1, ■ ■ ■ , 12) are some constants, whose expressions can be found 
in appendix E. By further using (3.34) and (3.36), we can reduce the above 
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expression as 



2 3 4 ^339) 



where Di, D2, Fi-F^ are some material constants and E is the Young's modulus. 
The formulas for these constants are given by 



_rfr7i + r|?7i + rir2(?7i - 8772) 
iJi — 



2r2(ri + r2) 



1 



or^ 1^1 ^2 J 



ri^2) - 4rfr2(-3r/2 + rie2) + 6r2r|(-r/ir/2 - 2r/2 + n^s)], 

-c/ — , -Ti — , r2 

F3 



Ti - r| n + r2 2ri + r2 
, ri — , r2 — • 

Ti 6ti 6ri 



12r3(ri + r2 

+ T^T2i2T2 - 77]l + 12r]2)]. 

By retaining the original dimensional variable, we obtain 



^ =2aE [W^ + ^D,V^ + ^D2V^ + a\F,V^ + F2VVxx] 



(3.40) 



+ a\F^VV^ + F,V^Vxx)l 



where V{X) = euo^ = Uqx is the leading order term of the axial strain. Then 
the total elastic potential energy of the layer should be given by 

^E= [ ^dX =2aE [ llv^ + Id^V^ + \d2V^ + a^{FiV], 

Jo Jo 2 6 4 (^3_42) 

+ F2VVXX) + a\F,VV^ + F^V'Vxx)]dX. 
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Although the expression oi given in (3.41) is one-dimensional, it is derived 
from a higher-dimensional setting and has the higher- dimensional effects built in. 
It should be noted that the expression (3.41) also contains higher-order derivative 
terms which play the role of regularization. Trunskinovsky (1982, 1985) (see 
also, Aifantis & Serrin 1983 and Triantafyllidis & Aifantis 1986) pioneered the 
idea of the regularization augmentation for solid-solid phase transitions which 
involves adding terms (like a strain gradient) into the usual constitutive stress- 
strain relation. Here the gradient terms in our expression is also derived and its 
coefficient is explicitly given. 

4 The mechanical dissipation functions 

One of the most important experimental results is that the measured en- 
gineering stress-strain response is rate-independent hysteretic. To model this 
rate-independent hysteresis phenomenon, one needs to take into account the me- 
chanical dissipation due to phase transitions. In this section, we shall derive the 
expressions of the mechanical dissipation functions in terms of the axial strain 
for the purely loading and purely unloading processes, which can be used to 
determine the total amount of energy dissipated during the phase transitions. 

The constitutive form of the mechanical dissipation rate ( has already been 
given in (2.6). From (2.6), it can be seen that the energy dissipation only occurs 
whenever a ^ and also depends on the phase state variable a. Thus, to 
determine the energy dissipation rate at certain material point, we first need to 
determine the value of the phase state variable. 

Based on assumption 1 and the phase transition criteria (2.7)-(2.8), we have 
already proposed the evolution laws of the phase state a for the purely loading 
and purely unloading processes in (2.9) and (2.10). 
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For the purpose of simplicity, the mechanical dissipation effect will be con- 
sidered in a one-dimensional setting in this paper, i.e., we neglect the influence 
of the radial strain on the mechanical dissipation. In this case, the deformation 
gradient tensor F should be replaced by the axial strain V. Thus, the evolution 
laws (2.9) and (2.10) reduce to 



a+{V) = { 1 



and 



a-{V) 



where the critical strains Vq^ and V-^ are determined by 



0, 


if 




< 1, 


if 




Oiy, 

V 


if 


<y < Vx 








0, 


if 




' 1, 


if 


Vx <y < y\ 




if 


<y < y{ 



>9^(K)+,«) 
da 

5^(K,-,«) 



da 



„=o = ^+(0), 
-A-(0), 



|q=0 



9^(V,a) 
da 
d^{y{,a 
da 



U=i = -^~(i) 



and ay satisfy 



(4.1) 



(4.2) 



(4.3) 



d^>{y,a) 
da 



A+(a+), 



di>{y,a) 
da 



-A {ay). 



(4.4) 



Here, y* denotes the maximum strain value in the experiments. It is easy to see 
that a^iy) are continuous functions. In principle, once ^, A'^{a) and A~{a) are 
given, we can obtain the explicit expressions of the phase state functions ay. 



In the paper of Rajagopal & Srinivasa (1999), some specific constitutive 
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structures were proposed for the Helmholtz free energy "^{V^a) and the dissi- 
pative response functions A'^{a) and A~{a). Through some calculations, it was 
found that a'^iV) were monotonically increasing functions during the intervals 
{Vf^,V-(^). Based on the model proposed by Rajagopal & Srinivasa (1999), we 
assume that a^{V) have the following property: 

Assumption 3. a'^lV) are monotonically increasing functions during the 
intervals {Vf^,V^). 

By virtue of the phase state functions a~^{V) and a~{V), we can express the 
total amount of mechanical dissipations during the purely loading and purely un- 
loading processes as functions of the axial strain V{Z). Suppose that the loading 
process starts (say, at time to) from the homogeneous configuration V{Z) = 
{0 < Z < I). Then the mechanical dissipation (say, at time t) through the loading 
process should be given by 

fl i*a i*t 

/ (dtdXdY 

Jo J~a J to 
pi pa pt 

= 111 A+{a)adtdXdY 

Jo J~a J to 

rl pa pa+{V{X)) 




pa' (V(A)) 

/ A+{a)dadYdX (4.5) 

Jo 



'0 J -a JO 



2a I I A+(a+(t;)) ' dvdX 



Jo 



2aE [ (l)j{V{X))dX, 
Jo 



where 



-{V) = ^£ A^a^iv))^^dv. (4.6) 



Similarly, for the unloading process (suppose that the unloading process 
starts from the homogeneous configuration V{X) = V* > Vf (0 < X < I)), 
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we can deduce 

$^ = 2aE f rAV{X))dX, (4.7) 

where 

^-AV) = -^^X"^^ ("C"^)^^^^- (4.8) 

The two functions 0^(V^) and 0^(V) defined in (4.6) and (4.8) are referred as the 
dissipation density functions. 

In principle, once A^{a) and A~{a) are given, (p'^iV) and <p~^{V) can be 
obtained. We also point out that since the strain is a measurable quantity, it may 
be easier to determine 0^(V^) and 0^(^) experimentally than yl+(a) and A~{a). 



5 The equilibrium equation and analytical solu- 
tions 

With the expressions of the total elastic potential energy and the mechanical 
dissipation functions, we can determine the equilibrium configurations of the 
layer during the phase transition process. We shall further consider an illustrative 
example with some given material constants and some special forms of dissipation 
functions in this section. Subject to the free end boundary conditions at the two 
ends of the layer, we construct the analytical solutions, which can capture some 
important experimental features. 

5.1 Equilibrium configurations 

In this subsection, we shall determine the equilibrium configurations of the 
layer during the phase transition process by using the variational method. 

In a general case, we can prove that for the purely loading and purely unload- 
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ing processes, the equilibrium configuration of tlie SMA body can by determined 
by using tlie principle of maximizing total energy dissipation (see Wang & Dai 
2009). The "total energy dissipation" is here referred as the part of work done by 
the external force that is not converted into the elastic energy or used to overcome 
the mechanical dissipation due to phase transformation. 

In the current case, the total energy dissipations during the purely loading 
and purely unloading processes should be given by 

=2aE I i{V{X) - 0)dX - {^e{V) - $^(0)) - ^+{V) 
Jo 

=2aE [\jV - il-V' + l-D,V' + ]d,V' (5-1) 
Jo 2 6 4 

+ a\F,Vl + F^VVxx + F^VV^ + F.V'Vxx)) - <Pj{V))dX. 

and 

<^-{V) =2aE [ ^{V{X) - V*)dX - {^e{V) - - ^^iV) 
Jo 

=<^*E-2aE I -fV*dX 

/° (5.2) 

+ 2aE ^ i^V - i^V^ + ^DiV + ^D^V^ 

+ a^FiV^ + F^VVxx + F^VV^ + F^V^Vxx)) - 0^ (^))^^, 

where 7 = P/ {2aE) is the dimensionless engineering stress acting on per unit 
reference area of the cross-section of the layer and represents the total elastic 
potential energy of the layer corresponding to the initial state of the unloading 
process. 

To determine the equilibrium configurations of the layer during the phase 
transition process, we need to find V{X) such that (S'~^(y) and (S'~{V) attain the 
maximum values. From (5.1) and (5.2), by using the variational principle, we 
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obtain that the equihbrium configurations should satisfy the following equation 
V + D^V^ + D2V'' + a'i-^Vxx + D^Vl + 2D^VVxx) + = 7, (5-3) 

where D-^ = 2F4 — F3. We refer (5.3) as the equilibrium equation. 

Remark: From another point of view, the equilibrium configurations of the layer 
during the phase transition process can also be determined by using the princi- 
ple of minimizing the total "pseudo-potential energy" of the whole mechanical 
system. The concept of "pseudo-elastic energy" function was first introduced by 
Oritz & Repetto (1999) to study the dislocation structures in plastically deformed 
crystals. 

In our case, the pseudo-elastic energy functions W^{V) corresponding to 
loading and unloading processes should be defined by 

W+{V) = ^e{V) + ^+{V), (5.4) 

and 

W-{V) = ^e{V) + $^(r). (5.5) 

Based on the pseudo-elastic energy functions, the total pseudo-potential en- 
ergies of the whole mechanical system can be defined by 

n^ = W^-P [ VdX. (5.6) 
Jo 

With the expressions of the total pseudo-potential energy (5.6) and by using 
the variational method, we can also derive the equilibrium equation (5.3). 

In fact, from (5.1), (5.2) and (5.6), it is easy to see that the principle of 
maximizing the "total energy dissipation" is equivalent to the principle of mini- 

29 



mizing the "total pseudo-potential" energy. For the purpose of convenience, we 
shall only use the principle of minimizing the total pseudo-potential energy in the 
sequel. Notice that both the principles of maximizing the total energy dissipation 
and minimizing the total pseudo-potential energy do not reveal anything about 
the path taken by the layer from its initial configuration to the final equilibrium 
configuration, i.e., these principles are entirely silent about the process and are 
only concerned about the initial and final states. 

Remark: Mielke et al. (2002) proposed a rate-independent, mesoscopic model 
for the hysteretic evolution of phase transformations in SMAs. In their model, 
an extremum principle, called the postulate of realizability (cf. Levitas 1995a, 
b), was adopted to determine the stable state of the material during the phase 
transformation process. We can also demonstrate that the extremum principles 
we used in this section are consistent with the postulate of realizability (cf. Wang 
& Dai 2009). 

5.2 The analytical solutions 

In this subsection, we shall construct the analytical solutions for an illustra- 
tive example with some given material constants Z^i, Z^2, -D3 and some special 
form of dissipation density functions (t>^{V). It will be seen that the solutions 
obtained can be used to explain some important experimental results. We also 
point out that most of the conclusions drawn from this illustrative example also 
hold in the general case. 

First, without loss of generality, we take the length of the layer to be 1. At 
the two ends of the layer, we impose the free end boundary conditions. By 'free 
ends', we mean that 

Vx = ^, at X = 0, 1, (5.7) 
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which is sometimes called natural boundary conditions and has been used by 
many authors (e.g., Ericksen 1975; Tong et al. 2001). 

Next, we shall give some further discussion on the dissipation density func- 
tions (p^iy). Based on (4.6) and (4.8) and by considering the properties of the 
phase state functions a^{V), we get that 0^(^) have the following general prop- 
erties: 

• (l)~^{y) is a continuous monotonically increasing function; (p^iV) is a con- 
tinuous monotonically decreasing function. 

• (t)'l{V) equals for < < and is a constant for <V <V*- (t)'2{V) 
equals for V{' <V <V* and is a constant for < V < Vq . 

Due to the diversity of materials, the dissipation density functions </'^(V^) 
could have many different forms. In this section, we choose 4>^{V) to be fourth- 
order polynomials in the intervals (Vj^, V^), which are given by 

^iiV) = HtiV - O + H^V - + H^V - V,^r + HtiV - V,^)' 
= + HtV + H^V^ + H^V^ + H^V^ 

(5.8) 

and 

0,-(V-) = H^iV - Vr) + H^{V - V,-f + H^{V - V,-f + Hl{V - V,-f 
= Ho + HiV + H^V^ + H^V^ + H^V\ 

(5.9) 

Remark: Here 4>^{y) are chosen to be fourth-order polynomials for the purpose 
of simplicity. In the general case, if 0^(V^) are smooth enough, we can also 
consider the Taylor series expansions of (t)^{V) in the intervals (Vq^, V^). 
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By substituting (5.8) and (5.9) into the equilibrium equation (5.3), we obtain 

-^ + V + D,V^ + D^V^ - ^a''Vxx+a''{DsV^ + 2DsVVxx) = 0, 

if < V < V^o^ or V^^ <V < V\ 

-(7 - ^0^) + H^V + HfV^ + HfV^-^a'Vxx + a\D,V^ + 2D3VVXX) = 0, 

if Vq^ <V < V{^, 

(5.10) 

where 



H^ = Hf, Hf = l + 2H^, H^ = Di + 3H^, Hf = D2 + AH^. (5.11) 

For the purpose of illustration, we choose the following numerical values in 
this paper 

Di = -23.81, D2 = 158.73, D3 = -20/3, V* = 0.1 

i^^o^ = 0.0087, ifi+ = 0.3132, = -7.1069, ^3+ = 41.9287, (5.12) 

Hq = 0.0097, H{ = 0.1497, = -5.4717, ^3" = 41.9287. 

With the above chosen material constants and through some calculations, we can 
get the values of some critical strains and stresses, which are given by 

Vq+ = 0.03, 1/+ = 0.083, Vq- = 0.017, = 0.07, 

^+ = 0.00973603, ^2^ = 0.01285714, ^+ = 0.01143011, (5.13) 

^- = 0.00777778, ^2" = 0.01089889, = 0.00920481, 

where ^+ are the valley stress values corresponding to V-{^, ^+ are the peak stress 
values corresponding to Vg^ and ^+ are the Maxwell stress values. 

Next, we shall solve equation (5.10) under the free end boundary conditions 
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(5.7). Here, we omit the detailed derivation process, which can be found in our 
previous work (see Wang & Dai 2009), and just give the expressions for the 
analytical solutions. 

First, we find that there exist constant solutions, which are given by the real 
roots of the following equations 

V + DiV^ + D2V'\ Q<V <V^ or <V < V\ 
7 = { (5.14) 
+ HfV + H^V^ + H^V^, <V < Vt- 

The constant solutions correspond to the homogeneous deformations of the layer. 
It's clear that these constant solutions satisfy Vx = for X G [0,1], thus the 
boundary conditions (5.7) are satisfied. It should be noted that as 7 varies, the 
number of the constant solutions can also be different. In fact, there is only one 
constant solution if < 7 < .^f or 7 > there are two constant solutions if 
7 = ^j*^ or 7 = S,2] there are three constant solutions if < 7 < S,2 ■ 

Besides the constant solutions, there are also non-trivial solutions when < 
7 < The expressions for the nontrivial solutions are given by: 




(5.15) 
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where is an integration constant and 



fHV) = -iV + ]^V^ + ^-D,V^ + -^D^V\ Q<V< V^, 

ftiV) = -iV + ]^V^ + ^-D.V' + -^D^V^ + M2±, < r < y*. 

(5.16) 



gf and (7^ are the two real roots of the equation 



+ fHv) = 0. 



(5.17) 



Notice that for non-trivial solutions to exist, equation (5.17) must have four real 
roots oti < gf < g2 ^ ctf. The constants and given in (5.16) are 



determined by M+ = H+, M+ = H+ + H+V+ + H^V+^ + H^V^^^ + ^4+ V^ 



Mf = Hq + M2" and = -{H^ + H^Vq' + H^Vq 



- ifg-Vo + ^rK) ) 

such that /i^(K)^) = /2^(K)^) and ft{Vi^) = /3^(^i^), which imply that Vx is 
continuous at V^"^ and V-{^. 

The constant is determined by the following equation 




n 



1,2,3,- •• . 



(5.18) 



Once is known, from the above relationship and through some calculations, 
the corresponding non-trivial solution can be obtained from (5.15). In this paper, 
we only consider the non-trivial solutions v^^ and corresponding to n = 1 and 
= 2 (it can be shown that the non-trivial solutions corresponding to large n 
cannot be the preferred solutions). 

We choose the half thickness-length ratio a = 0.00866. For the purely loading 
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Figure 1: The 7 — A curves for the constant solutions and the first two non-trivial 
solutions in the loading process. 

and purely unloading processes, we plot the stress-strain (7 — A) curves corre- 
sponding to the constant solutions and the first non-trivial solutions in figure 1 
and 2, respectively. In figure 1 and 2, the three dashed curves labeled vf, vf and 
vf correspond to the three constant solutions and the two solid curves labeled 
v^i and v^2 correspond to the two non-trivial solutions. 

From figure 1 and 2, we can see that when A^ < A < Af, there exist multiple 
solutions. To determine which solution is the preferred one, we compare the 
total pseudo-potential energies for all the possible solutions. In the displacement- 
controlled problem, the total pseudo-potential energies for the purely loading and 
purely unloading processes should be given by 



--2aE 



3 ^ 4 



T:a'VVxx 
b 



(5.19) 
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Figure 2: The 7 — A curves for the constant solutions and the first two non-trivial 
solutions in the unloading process. 

and 

^- =w-(v) = ^E{y) - + ^{v) 

=2aE [ il-V' + Id,V' + ^D^V - ^a^VYxx (5-20) 
Jo 2 S 4 6 

-a'D,VV^ + <f)-{V))dX-<^%. 

We plot the differences of the pseudo-potential energy Al]^ between the first two 
non-trivial solutions and the constant solutions in figure 3 and 4. From figure 
3 and 4, we can see that for A < A^ or A > A^, the constant solution is the 
preferred solution and for Af < A < A^ the first non-trivial solution v^^ is the 
preferred one. 

In figure 5, we show the 7— A curves corresponding to the preferred solutions, 
where the curve labeled represents the loading part and the curve labeled 
"— " represents the unloading part. We can see that the curves shown in figure 
5 capture some important features of the engineering 7 — A curves measured in 
experiments (see Shaw & Kyriakides 1995, 1997; Sun et al. 2000; Tse & Sun 
2000). For example, the stress peak (for the loading part) and stress valley (for 
the unloading part), the stress plateaus, the hysteresis loop and so on. 
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AQ^ (2aE) 




Figure 3: Differences of the pseudo-potential energy between the non-trivial 
lutions and the constant solutions for the loading process. 



AQj (2aE) 




Figure 4: Differences of the pseudo-potential energy between the non-trivial 
lutions and the constant solutions for the unloading process. 
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Figure 5: 7 — A curves for the preferred solutions. 



To draw figure 5, we just use the principle of minimizing the total pseudo- 
potential energy. We suppose that the configuration of the layer can jump from 
one metastable solution to another metastable/stable solution once the pseudo- 
potential energy of the latter one becomes smaller than the former one. But 
in fact, there may exist potential barriers between two metastable solutions. In 
this case, if we don't give the layer enough perturbation, the layer will keep 
in its original metastable solution until this solution become unstable or some 
other limitation conditions happen. Thus, the curves shown in figure 5 are not 
necessarily the actual curves measured in the experiments, as metastable solutions 
may be maintained. 

Alternatively, we can use the "limit-point" instability criterion for the onset 
of the nucleation or coalescence process. For the loading part, we assume that the 
layer keeps in the first constant solution until the total displacement A reaches 
Vq (a "limit-point"; see figure 1), which corresponds to the peak stress value. As 
the total displacement goes on increasing, the martensite phase will nucleate even 
for the homogeneous deformation, thus the configuration of the layer jumps from 
the constant solution to the first non-trivial solution, i.e., the nucleation process 
happens. As A further increases, the deformation follows the first non-trivial 
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Figure 6: 7 — A curves based on the "limit-point" criterion. 

solution. Once A reaches A^, another "limit-point" is also reached (see figure 1), 
the deformation of the layer jumps to the third constant solution. Thus, it is at 
the displacement value A = A4 that coalescence process happens. Similarly, for 
the unloading part, we assume that the nucleation process happens at A = V^i (a 
"limit-point") and the coalescence process happens at A = A3 (another "limit- 
point"; see figure 2). The corresponding 7 — A curves are shown in figure 6. We 
can see that the curves shown in figure 6 are consistent with the experimental 
results (see Shaw & Kyriakides 1995, 1997; Sun et al. 2000; Tse & Sun 2000). 

6 Analysis on the coalescence process 

In this section, we shall give some analysis on the coalescence process, which 
plays a central role in the whole phase transition process. The following analysis 
of this section are based on the analytical solutions obtained in the section 5. 
Here we wish to give some descriptions and explanations for the origin of the 
instability during the coalescence process, the accompanying stress drop/jump 
and the morphology varies of the specimen. We shall also consider the size-effect 
of the specimens on the coalescence process. 
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As mentioned before, the coalescence process of phase fronts observed in 
quasi-static experiments is a dynamic one, which is an indication that an in- 
stabihty occurs. In subsection 5.2, we have aheady used the minimum energy 
principle to model the coalescence process (see figure 5). By comparing the total 
pseudo-potential energies of the constant solutions and nontrivial solutions, we 
found that the coalescence process should take place at the point A = for the 
loading case and A = for the unloading case. Besides the minimum energy 
principle, another possible reason for the onset of the coalescence process could 
be the fact that some "limit-points" have been reached. In subsection 5.2, we 
also used the "limit-points" instability criterion to model the coalescence process 
(see figure 6). We found that it was at the points A = for the loading case 
and A = A^ for the unloading case that some "limit-points" had been reached 
and the coalescence process took place. 

Remark: In the sequel, we shall only use the "limit-points" instability criterion 
to model the coalescence process. 

Systematic experimental results have obviously shown that the coalescence 
process is inevitably accompanied the varies of the stress value and the surface 
morphology of the specimen. Here we shall use our model to describe these 
phenomena. For the purpose of clearness, we consider the second nontrivial 
solution (i.e., we choose n = 2 in equation (5.18)) and consider the case that the 
phase fronts coalesce at the middle part of the layer. Thus the transformation 
scheme can be simply written as MA + AM — >■ M for the loading case and 
AM + MA A for the unloading case. Figure 7(a) shows the engineering stress- 
strain curve corresponding to the coalescence process during the loading case 
(a = 0.0043). Corresponding to the 5 points A-E show in figure 7(a), we plot the 
profiles of the layer in figure 7(b). Here the radial deformation has been enlarged 
for clearness. From figure 7, we can see that the stress drop and the deformation 

40 



y 




Figure 7: (a) Engineering stress-strain curve for the coalescence process during 
the loading case (a = 0.0043). (b) Profiles of the layer corresponding to the 5 
points A-E. 
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Figure 8: (a) Engineering stress-strain curve for the coalescence process during 
the unloading case (a = 0.0043). (b) Profiles of the layer corresponding to the 5 
points A-E. 

process of the layer are very similar to the actual experimental procedures. Figure 
8 shows the stress jump and the deformation process of the layer corresponding to 
the coalescence process during the unloading case (a = 0.0043). We can see that 
the coalescence process shown in figure 8 is also consistent with the experimental 
results. 

Next, we consider the influence of the geometric size effect of the layer on the 
coalescence process. From (5.15) and (5.18), we can see the nontrivial solutions 
obtained here depend directly on the half thickness-length ratio of the layer. 
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Figure 9: 7 — A curves of the second nontrivial solutions corresponding to some 
different half thickness-length ratio (a = 0.0026, 0.0035, 0.0043). 

which shows the fact that our model can reflect some important information on 
the geometric size effect. Here we only consider the loading case, and the similar 
result can also be derived for the unloading case. Figure 9 shows the 7 — A curves 
of the second nontrivial solutions [n = 2) corresponding to some different half 
thickness-length ratio. From figure 9, we can see that as a decreases, the 7 — A 
curve of the non-trivial solution moves towards the constant solution. Actually, 
one can prove that the 7 — A curve of the non-trivial solution moves toward first 
constant solution when < 1 < ^2 ^^'^ towards the third constant solution 
when C,i < 7 < ^m- The following is the proof for the case of the second nontrivial 
solutions {n = 2) and ■Ci^ < 7 < 'Cm (^^e case of other nontrivial solutions and 
< 7 < .^^ can be proved similarly). 
First, for the second nontrivial solution, we have the following equation (cf. 



which can be considered as a relationship between the half thickness-length ratio 
a and the constant C~^. Through some simple analysis, we find that the equation 



(5.18)) 
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+ f~^{V) = has four real roots < gf < g2 ^ (^t only if 

C2+<C+<min(C+,C3+), (6.2) 

where 

Ct = -fHvt), ^ = 1,2,3. 
In this case, the function + f~^{V) can be written in the following form 

+ f^{V) = fc{V){V - at){V - gt){V - gt){V - at), (6.3) 

where fc{V) is a non-zero and bounded continuous function depending on the 
constant C. 

Through some further analysis, we find that when C,i < ^ < C'^ should 
satisfy < C+ < and 

lim af < vf < lim g^ < < lim gt = '^t ~ ^ (6-4) 

c+^c+ c+^c+ c+^c+ c+^c+ 

where v^, and are the three real roots of equation (5.14). From (6.1) and 
(6.3), we have 



^ (6.5) 



9 / I/Q-D3V 

J at V fc{V)(V-a+)(V~g+){V~g+)(V~a+)'^'' 



Thus as a tends to zero, the integration f^? .1 — — +— - — T^dV 

-'ai \l fc{V){V-aJ){V-gJ){V-g:^){V-a^) 

should tends to infinity, and this is equivalent to C+ tending to . By using 
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(5.15) and (6.5), we can give the following derivation 



VdZ = 2 / VdZ = 2 / V—dV 

Jo Jg+ dV 
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fciV)iV-at)iV-gi)iV-gt)iV-at)'^^ 

raty I I/G-D3V ,y 

•'gt \J fc{v)(v-o^t)iy~9tw~9t){y-^t) 



•'gt V fc(V)iV-a+)iV-g+)(V-g+)(V-a+r'' (O.b) 
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Jgf ^ W fcivw-»t){v-9t)iv-gt)iv-4r^ 

rgt I l/a-DzV jy 

-fgi V fc{V){V-a+){V-g+)iV-g+){V-a+r'^ 



[92 / il/6-D3V){g+-V) ,y 

V Y /cW{V-a+)(V-g+)K-^) 

[92 I I/6-D3V jy ' 

Jn+ \ fr,(V\(V-n+MV-a+MV-n+\(V-n+-\ 



'gt V fc{V){V-a+){V-gtW-g+){V-a+) 

By using (6.4), it is easy to prove 



r il/6-D3V)igt-V) 



and 
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y I I 1/6 - D3V , 

c%Jgt y fc{V){V-at){V-gt){V-gt){V-atr 



Thus, as a tends to zero, the total elongation of the layer should satisfy 



limA= lim A= lim (7^ = f-^- (6-7) 



c->c+ c->c^ 



This means that the 7 — A curve of the second nontrivial solution tends to that 
of the third constant solution as a tends to zero when ■Ci' < 7 < 'Cm- 

Figure 9 shows that the 7 — A curve of the nontrivial solution has a snap- 
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I A 

Figure 10: Engineering stress-strain curves for the coalescence process during the 
loading case (a = 0.0026, 0.0035, 0.0043). 

back structure. But in a purely loading process, the total elongation A cannot 
decrease. Thus the state of the layer should jump from the nontrivial solution 
to the constant solution at some special point, which corresponds to the coales- 
cence process. Here we also use the "limit-points" instability criterion to model 
the coalescence process. The actual 7 — A curves of the coalescence processes 
corresponding to some different half thickness-length ratio are shown in figure 10. 
From figure 10, we can see that as a decreases, the stress drop during the coales- 
cence process also decreases. Actually, (6.7) implies that the limit point should 
tend to the third constant solution as a tends to zero. Thus, the stress drop 
should tend to zero as a tends to zero. Similarly, in the process of AM+MA^A 
during the unloading case, the stress jump should tend to zero as a tends to zero. 
Thus we can say that the geometrical size of the specimen plays an important 
role in the coalescence process. 
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7 Conclusions 



In this paper, we aim to study the phase fronts coalescence phenomena during 
the phase transition process in a thin SMA layer. For that purpose, we derived a 
quasi-2D model with a non-convex effective strain energy function while taking 
into account the rate- independent dissipation effect. 

In a two-dimensional setting, we studied the symmetric deformation of a thin 
SMA layer. Based on the field equations and the traction-free boundary condi- 
tions, by using the coupled series-asymptotic expansion method, we expressed the 
total elastic potential energy of the layer $£; as a function of the leading order 
term of the axial strain V{X). We further considered the mechanical dissipation 
effect in a purely one-dimensional setting and expressed the total amount of me- 
chanical dissipations $J (loading case) and $^ (unloading case) as functions of 
the axial strain V{X). The equilibrium equation was then determined by using 
the principle of maximizing the total energy dissipation. We further considered 
an illustrative example with some given material constants and some special form 
of dissipation density functions. With the free end boundary conditions, we man- 
aged to construct the analytical solutions for both a force-controlled problem and 
a displacement-controlled problem. 

Based on the analytical solutions obtained and by using the limit-point in- 
stability criterion, we studied the phase fronts coalescence phenomena. It was re- 
vealed that during the coalescence process, the configurations of the layer switched 
from the nontrivial solution modes to the trivial solution modes, which was caused 
by the presence of the "limit points" . The morphology varies of the layer and the 
accompanying stress drop/jump during the coalescence process can be described. 
The influence of the thickness-length ratio of the specimen on the coalescence 
process was also studied. It was found that the zero limit of the thickness-length 
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ratio can lead to the smooth switch of nontrivial modes to trivial modes with no 
stress drop or stress jump. 

Appendix A Incremental elastic moduli 

For initially isotropic material, in the case that there are no prestresses, the 
elastic potential energy function should be a function of the principal stretches 
Ai and A2, namely ^=^(Ai,A2). Denote by = ||^|ai=A2=i, then ^1 = ^2 
should vanish since there are no prestresses. 

The non-zero first order incremental elastic moduli can be written as 

n = ai2i2 = 2^^^~ ^2), 
T4 = a\22i = Ts- 

There are only two independent constants among Tj. 

The non-zero second order incremental elastic moduli can be written as 





2 

~ '^111111 


= ^111, 




V2 


2 

~ '^111122 


= ^112, 




V3 


2 

~ '^111212 


1, 

= -^{Ti + T2 + r]i - 


V2), 


V4 


2 

~ '^111221 


1, 

= -^{-n -T2 + r]i 


-V2), 



There are only two additional independent constants among rji. 
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The non-zero third order incremental elastic moduli can be written as 



6'i = 

^1 


"11111111 


^ 1111 1 


02 = 


illillzz 


^1112, 


03 = 


11111212 — 


^(-3ri - 3t2 + 3r/i + 3r/2 + 29i - 262), 


04 = 


'-^1 1 1 1 1 991 — 


-^(3ri + 3r2 - 3r]i - 3r]2 + 29^ - 262), 


05 = 


'-^1 11 19999 ~ 


^11225 


a 

C76 — 


- 

"11121222 ~ 


— dri — 6T2 + D?72 + ^1 + /:y2 — Ot/5), 


07 = 


"^^11122122 = 


^(3ri + 3r2 - 67/2 + ^1 + 2^2 - 3^5), 


0s = 


"12121212 ~ 


+ 3x2 + 6r/i - 6r/2 + ^i - 4^2 + 3^5) 




09 = 


"12121221 ~ 


^(-3ri- 3x2 + ^1-4^2 + 3^5), 



010 = 


- "12122121 - 


= ^(3ri + 3r2 - 2r]i + 2r]2 + 0i - 4^2 + 3^5 



There are only three additional independent constants among 6i. 
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Appendix B Non-dimensional field equations 

The full forms of the non-dimensional field equations (3.15) and (3.16) are 
given by: 

2t3Us + (r2 + T3)w.j; + s{Ar3Uss + (2r2 + 2t3)w^^s) + nu.^^ + t{2ri^wus + 2r]3UsU^ 

+ (4?72 + 4:r]i)wsW^s) + V2wu^x + Vi^xUxx + s{16r]3UsWs + ^^^wu^s + Ar]3UssU^ 
+ {2r]2 + Sr]i)wsW^ + 8773^^^^^ + 4?74ii;^.u^.s + (2?72 + 2r]^)ww^s + (2?72 + '^Vi)uxW^s 



+27]2WsU^x + 2r]^UsW^^ + r]3WxW^^)) + e^^Osiu'^Us + 29gwUsU^ + 6'3nsM^ 

6 6 6 6 

+ (y + ^)w^^w^x + (6'5 + 67)wu^w^ + S^(86'3W^'«ss 

+ 166'3M^W^M;ss + WiWsWssWx + (46^2 + 46'4)w^Wa;s) + ^^'sW^M^^ + 62WUxUxx 

+ + S^(286'3M^W^ + W^ulUss + W'iWWsUss + W^WUsWss + S6qWsUssU^ 

+ 86qUsWssUx + (26*2 + 146*4)^^11;^ + 86gUsUssWx + 46'4WWssWa; + A6iWssUxWx 

+ 26'ion^sW^^ + 166'6UsU;sU^^ + WjWsWxUxs + (46^6 + 46'9)u^ty3;3 + (46*2 + 46'4)u;w;stya:< 
+ (46*5 + A6j)wsUxWxs + (46110 + 4:6j)usWxWxs + {6q + 6q)wIwxs + 26^w'luxx 

'U^x'U^xx 

) + S(46'8M^ + 166'3WM3W^ + 26310^ Uss + 166'6MsW5Mx 
+ 4^6WMssMx + 263UssUl + (2^6 + e^g)^^^^ + (2^2 + 8^4)^^,^^ + (3^io + 267)UsW 

6 6 

+ (26^5 + 867)WsUxWx + + + ^O&WUsUxs + 863UsUxUxs + 4:67WWxUxs 



+ A6iUxWxUxs + (6*2 + 6'4)m;^w^s + (26^5 + 267)wUxWxs + {62 + 6i)ulwxs + 26'3M^M; 



2., 



+ 26^WWsUxx + 262WsUxUxx + 26aUsWxUxx + ^63wlUxx + 26iWUsWxx + 26aUsUxWx: 
+ 6qWWxWxx + 63UxWxWxx)) = 0, 

(51) 
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GtiWs + {2t2 + 2t3)u^s + TsWxx + 4:STiWss + e(4?73M^ + GrjiwiUs 

+8s^r]iWsWss + Qr]2WsUx + Qri^UsW^ + 2?73tf;i^. + (2?72 + 2rii)wUxs 

+ {2r)2 + 2'r]4)u^Ur,s + 2r]4UsU^x + VsW^Uccx + Vs^^xx + VsUxWccx 

+s(12?7im;^ + SijsUsUss + 4:r]iWWss + "^^WssU^ + Ar]iUssW^ 

+ (4?72 + 4?74)ty,u^^ + 8r]iUsW^s + 'irjsyJxWxs + 2r]3WsW.^.^)) 

+ e^{4:93wu1 + SOiw'^Ws + Ss^Oiw^Wss + 4:dQulux + 692WWsUx + 36'5ty, 

+ 66'4WMsii'2; + 69jUsUxWx + 293wwl + 29qUxwI + (^2 + 9i)w'^Uxs 
+ (26^5 + 297)wUxUxs + (6'2 + 9i)ul.Uxs + 297WUsUxx + 29iUsUxUxx 

1 1 

+9qWWxUxx + 93UxWxUxx + -9'iw'^Wxx + 9qWUxWxx + -9zu'iWxx 

+ s'^{129iwl + IG^sUslfsU^s + 86'3U^tyss + WiWWsWss + 892WsWssUx 

+ 89iWsUssWx + WiUsWssWx + 293WssWl + (46'2 + A9A)w'^gUxs 

+ 169aUsWsWxs + W-iWsWxWxs + 293wlwxx) + s(206'3M^Ws + 126'iWW^ 

+ 86'3W;MsMss + 29iw'^Wss + l292wlUx + WqUsU^sUx + ^92WWssUx 

+ 295WssUl + 2A9aUsWsWx + A9aWUssWx + A97UssUxWx + 793WsWl 

+ {496 + 4^9)w^Mxs + (4^2 + 4^4)WW,M^, + (4^5 + 4^7)m;,m^m^,+ 

(46'io + 4:97)UsWxUxs + (6*6 + 9g)wlUxs + W^WUsWxs + 897UsUxWxs 

+4:93WWxWxs + 4:9qUxWxWxs + A97UsWsUxx + 296WsWxUxx + 29i0ulwxx 

+ 293WWsWxx + 296WsUxWxx + 29qUsWxWxx + lr98wlWxx)) = 0. 
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Appendix C 



The formulas of the constants [i = I,-- - ,25) in (3.26)-(3.28) are given 
below: 



ai 



a2 = - 



an 



ai6 



112 r2Jh Vs Va. 

'2ts 2ts^ 2x3 2r3^ 

V2 TjJh V3 V4 

2t3 2r3^ 2t3 2n 



V2 , n?73 
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2t, 2x3^' 



= — 



2t3 2x3^' 



^ mm _ T2m'^ _ mm 
~ 2x32 2x33 2x32 2x32 

V2V3 r2m^ vs^ , 'n3V'i 



T-32 r32 

V2V3 r2m'^ , V3V4 



2x32 2x3^ 2x32 + 2x32 

_ mm _ nm^ 'nh _ 0^ 

""'"2x32 2x3^+4x32 4r3' 

_ mm mm _ n^s^ _ ^ + Ii^ 

" 2x32 + 2r32 r,^ 2r, + 2x32 

_ mm nv3^ , nOs 

O'lo — 7i — n — — n — r 



0^ 


To 6*3 


_ A 


0A 


4r3 


4- 4 
4r32 


4r3 


4X3' 




2x32 + 


Oe 


^7 


2r3 


2r3 " 


"2r3' 




r,^3 






4r3 


4r32 


4r3 


4r3' 



2x32 2r3^ 4r3 4r32 ' 
V2_ ^ r2m _ ^ _ T-2^?73 _ 2r2r/3 _j_ ^ 2r2?74 
6ri 6riT3 Sn Snr^^ Sms Sn StiT-s' 



a - ^2 I '^2^2 ^3 T2773 r2^773 773 r2?73 ^ m ^ Vi ^ T2m 
3ri 3rir3 3ri 2x32 6x1X32 2x3 3^X3 3ri 2x3 3^X3^ 



^1 _^ r2r/i _ ri?73 _ T2?73 _ V3_ ^ m_ 
6ri 6rir3 6x32 6x32 6x3 6x3 ' 



3ri2 6ri2r3 3ri 3rir3 3ri 6x1X32 3x1X3 3xi 3x1X3' 

^ _ m__ TjJh ^ T2m _ 'T2^m _ V3_ _ r2^m _ ^2^3 _^ ^ _^ ^2^74 



3xi 3xi2 3x1X3 6xi2x3 3xi 6x1X32 3x1X3 ' 3xi ' 3x1X3' 

m _ T2m J]2_ ^2^2 _ 72^73 m ^ Vi 

6x1 6x1x3 6x1 6x3 6x1x3 6x32 6x3 6x3 



^ _ vi ^ T2m m I m ^2^2 T-2773 773 ^ 774 
6x1 6x1x3 6x1 6x3 6x1x3 6x32 6x3 6x3' 
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12ri2 ^ 24ri2r32 ' SnVg 6ri ' 8x3^ ' 24r32 24^x32 ' 24r: 



3 

5r2?72 , 5?73 rir2?73 5x2^7/3 r2^?73 ri?73 5r2?73 , r2^?73 



7?73 _^ 13r2?73 _ 5?74 _^ ri?74 ^ 7r2?74 13r2^?74 ^ 7?74 13r2?74 



24x3 12rir3 12ri 8x3^ 24x3^ 24x1X3^ 24r3 Unr^ 

_ T2771 771 _ V2_ T2^2 Ti!/2 _ 5r2^772 T2^r]2 _ 5r2?72 

"24x3^ 24r3 6ri 12ri2 ^ 8x3^ 24x1X3^ 24ri2r32 12rir3 
T2^'72 , V3 n'^V3 7rir2?73 , r2^?73 , r2^?73 7rir]3 r2773 



+ 



6ri 24r33 24x3^ 4r33 12x1X33 24x32 2x32 

^2^'73 , ^73 , ^2% V4 , ^2^74 T-2^'74 Vl ^2^74 



O20 



4X1X32 8x3 3X1X3 6x1 6x32 8x32 6X1X32 8x3 3X1X3 

Vi , T2^Vi , ^2?7i ?72 , Tir]2 T2II2 T2^r]2 r]2 r2r]2 



24xi 24x1X32 12x1x3 8x1 6x32 12x32 8x1X32 12x3 4x1X3 

J73_ _ n^V3 _ nT-2?73 7X2 ^?73 ?"2^^3 _ ijTTs 7X2?73 r2^?73 

12xi 12x33 4x33 24x33 24x1X33 4x32 12x32 8x1X32 

_^ 7773 _^ r2?73 774 , rir74 7x2?74 X2^?74 7?74 X2?74 



24x3 6x1X3 12xi 4x32 24x32 12x1X32 24x3 6x1X3' 

_ _ V]_: ^iVi _ ^2^^1 _ ^2^1 V2 _ r2r]2 T2^ri2 _ r]2 T2r]2 

8x1 6x32 8x1x32 4xiX3 24x1 12x32 24x1x32 12x3 12x1x3 
_ ^2773 5xiX2?73 r2^?73 5xi?73 72?73 , ^ _ ^;2?74 _ V4_ 

3x33 24x33 8x33 24x32 4x32 + 8x3 8x32 8x3' 

_ ^2^71 , ^2^?7l , ^2^^7l V2 , 1"1'72 , T2V2 ^2^12 , ^72 

O22 —77; n + TTa o o + 7, o TT; r 



12xi2 24xi2x32 8xi2x3 12xi 24x32 24x32 8x1X32 24X3 

T"2'72 , ^73 T-iX2?73 X2^?73 Xi?73 X2^?73 X2?73 r]i 
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4x1X3 12x1 12x33 12x1x33 12x32 4x1x32 4x1x3 12x1 

^ riT]4 T"2^?74 T"2?74 



24x32 8x1x32 4x1 X3 ■ 



1"2'7l , Vl V2 , T2r]2 , Xi?72 X2^?72 , X2^?72 r2?72 



023 =77:; — + 



24x32 24x3 12ri 12xi2 24x32 8x1X32 24xi2x32 4x1X3 
_^ T2^^2 _^ ^3 _ T1X2773 _^ r2^?73 xi?73 , X2^?73 , r2?73 774 
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Xi2x3 12x1 12x33 12x1X33 12x32 4x1X32 4X1X3 12X1 

Xi?74 _ X2^?74 _ T-2?74 

24x32 8X1X32 4x1x3' 

Vi . '^2^Vi . ^2^71 m . Tiri2 X2?72 X2^?72 772 



24xi 24x1X32 12x1X3 24xi 12x32 12x32 24x1X32 12x3 
r2?72 ri2?73 ^ X2^r73 ^ r^rjs _^ T2V4. V-i 



12xiX3 12x33 12x33 6x32 12x3 12x32 12x3' 

025 = - 



Vl , X2^?7l X2?7l 7/2 X2?72 , 'r2^?72 ?72 



24xi 12x32 24x1X32 12x1x3 24xi 12x32 24x1X32 12x3 

+ 



r2V2 _ ri^V3 7-2^773 T2773 ^3 _ 7'2^4 _ _rh_ 

12X1X3 12x33 12x33 6x32 12x3 12x32 12x3 
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Appendix D 



The formulas of the constants 6j [i = 1, ■ ■ ■ , 22) in (3.29)-(3.31) are given 
below: 



h 

h 

h 

h 
h 

b9 
bio 
bii 
bi2 
bi3 
bi4 
bi5 
bi6 



-Sri^ + T2\-3r]i + r]2) - rir2(3r/i + 2r/2) + n\-5T2 + Tr/a) 

6ri(ri-r2) 

5ri^ + 2r2^?7i + n\5T2 - 7r]2) + rir2(4r/i + r/2) 
6ri(ri-r2) 

5ri3 + ri^(5r2 + 2r/i - 97/2) + 2r2^r/2 + rir2(4r/i + r/2) 
6ri(ri-r2) 

-Sri^ + rir2(-6r/i + 7/2) + r2^(-3r/i + r/2) + ri^(-5r2 + 37]^ + 4r/2) 

6ri(ri-r2) 

+ ri^(r2 - 37/2) + T2'^i-Vi + ^2) + rir2(?7i + 2772) 
6ri(ri-r2) 

-n"'^ - 3rir2?72 + T2^{-y]i + 772) + ri^(-r2 + 771 + 2772) 
6ri(ri - T2) 

+ 7'i^(r2 + 7/1 - 4772) + 3rir27/2 + T2^{-r]i + 7/2) 
6ri(ri-r2) 

+ ^^(^2 - 37/1) + r2^(7/i - 7/2) + rir2(27/i + 7/2) 
6ri(ri-r2) 

Tl^ + T1T2 - Tir]2 + r27/2 





2ri 


-2r2 




+ T-iT-2 - 


-rir]2 + T2ri2 




n 


- T2 


_v 


+ 7-1^-2 - 


-Tir]i + T2771 




2ri 


-2r2 




+ riT2 - 


- ^1'72 + T-27/2 




2ri 


-2r2 




+ TlT-2 ■ 


-riri2 + T2r]2 




2ri 


-2T2 




+ T-ir2 - 


- ri7/2 + r27/2 




2ri 


-2r2 




+ rir2 - 


- n?7i + 7-2 r/i 



2ri - 2r2 

2ri^ + rir27/i + 2ri^(r2 - 7/2) + r2^(27/i - 772) 
6ri(ri-r2) 
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W + ri\AT2 + 5r/2) + nT2{Ar]i - 2r/2) + r2^(r/i + 772) 

6ri(ri - T2) 
2ri^ + T2^r]i + 2ri^(r2 - r/2) + rir2(2r/i - r/2) 
6ri(ri-r2) 

+ ri^(r2 - 3r/2) + r2^(-r/i + r/2) + nT-2fa + 2r/2) 
6ri(ri-r2) 

-ri^ - 3rir2?72 + r2^(-?7i + 772) + Ti^{-T2 + r]i + 2772) 
6ri(ri - r2) 

+ ri^(r2 + 771 - 4772) + 3rir2772 + r2^(-77i + 772) 
6ri(ri-r2) 

+ ri^(r2 - 37/1) + r2^(7/i - 7/2) + rir2(27/i + 7/2) 
6ri(ri-r2) 
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Appendix E 



Cl 



The formulas of the constants q (z = 1, ■ ■ ■ , 12) in (3.38) are given below: 
rf + riT2{T2 - 2r]2) + (2r2 - r/a) + r|(2r7i + 772) 



6(ri - T2) 



2 



^2 =7n ^^(2rl^ " 2rj^(^2 + ^2) + rl{-2r]i + 772) - nr^ivi + 3r/2) 

6ri(ri - r2)^ 

+ TlT2i-AT2+Vl+QV2)), 

2rf + 2rir2r7i + r2(2r2 - 3r/2) + r|^2 



C3 



3(ri - r2) 



2 



C4 =^;r^^— -^(2ri^ + ri'(-4r2 + r/i - m) - ri{r]i + r/2) 

+ 3rir2(-2r2 + r/i + 772) + TiTl{-7r]i + 3772)), 
C5 = g^^^ ^^^^a (-3ri^ + r2^7/i - ri2(2r2 + 771 - 67/2) + rir2(r2 - 2(7/1 + 27/2))) 

C6 = ^ , ^ -^(-2ri^(r2 - 771) - T2^rii + riV2(-2r2 - 27/i + 7/2) 

6ri(ri - r2)2 

+ riT2^{-rii + 772)), 
2ri^ + 2rir2 + 2r27/i - 3ri7/2 + r27/2 



C7 

C8 



-12ri + 12r2 
^2^71 - riri2 - r2?72 



6ri 



cg =T7rV- ^ (2^1^ + ^1^(2^2 -Til- 4r/2) + 2r2^(r/i - r/2) + rir2(?7i + AT2)) 

12ri(ri - r2) ' 

2ri2 + ri(2r2 + r/i - 4772) + r2(r/i + 2r/2) 
Cio = 



Cll 



12 (-n + r2) 
2ri7/i + r27/i - 3ri7/2 - r27/2 
6ri 



C12 =T7r^ 7 (2n^ + ri2(2r2 - 47/i - 7/2) + 2r2^(7/i - 7/2) + rir2(47/i + 7/2)) . 

12ri(ri - r2) ' 
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